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5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a vehicle motion model 
generating device for generating a vehicle motion model which 
represents amotion state of a vehicle and amethod for generating 
10 a vehicle motion model , and more specif ically to a vehicle motion 
model using a recurrent neural network. 

2 . Description of the Related Art 

Researches and developments of various techniques have 
been hitherto carried out in order to enhance the operational 

15 stability of vehicles. As one of these techniques is known 
a vehicle motion model generated by modeling the behavior of 
a vehicle based on the kinematic theory of a vehicle . The vehicle 
motion model is generated through experiments or simulations 
by modeling a yawing motion, a horizontal motion, a rolling 

20 motion, etc. of a vehicle under some running condition, and 
set based on equations for motion of the vehicle. The 
operational stability of the vehicle can be estimated by 
analyzing the behavior of the vehicle , that is , the motion state 
of the vehicle or a vehicle motion state based on the vehicle 

25 motion model. 
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Furthermore, a road surface friction coefficient is used 
for the operation of control parameters in order to perform 
vehicle control such as traction control , braking force control , 
torque distributing control or the like. Thus, there has been 
5 proposed a technique of using the vehicle motion model as 
described above and estimating the road surface friction 
coefficient based on the vehicle motion model and the actual 
vehicle motion state. For example, JP-A-8-2274 discloses a 
method of estimating the road surface friction coefficient by 

10 using adaptive control . In addition, JP-A- 10-242030 discloses 
a method of estimating the road surface friction coefficient 
by comparing a vehicle slipping angle estimated by an observer 
with a reference value on a high jx-road based on the vehicle 
motion model and a reference value on a low ^i-road based on 

15 the vehicle motion model. 

JP-A-4- 138970 discloses a method of estimating the motion 
state of the vehicle by using a neural network. Specifically, 
a side slipping angle and a yaw rate are estimated on the basis 
of a neural network in which an easily-measurable parameters 

20 (for example, a longitudinal acceleration, a lateral 
acceleration, a vertical acceleration, a stirring torque, a 
front -wheel rudder angle, a vehicle speed, a rear-wheel rudder 
angle, etc. ) in the vehicle parameters representing the motion 
state of the vehicle are set as inputs and difficultly-measurable 

25 parameters (a side slipping angle and a yaw rate) are set as 
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outputs. Furthermore, JP-A-6-286630 discloses a method of 
estimating the road surface friction coefficient based on the 
running condition of the vehicle detected by using the neural 
network. According to the JP-A-4-138970 and JP-A-6-286630, 
5 in the neural network, adjustment (learning) of the coupling 
weight coefficient is carried out in advance according to an 
algorithm such as back propagation so that the output corresponds 
to a teaching signal. 

According to the method described in the JP-A-8-2274 or 

10 JP-A- 10-242030 , when a vehicle motion model is set, a motion 
equation is linearly approximated to avoid cumbersome operation 
processing in the solution calculating process. Therefore, 
the vehicle motion model may not accurately represent the motion 
state of the vehicle, that is, the behavior of the vehicle in 

15 a non-linear region. 

Furthermore, according to the method described in the 
JP-A-4-138970, a feed-forward type neural network is used, the 
value output from the neural network and the value input to 
the neural network are independent of each other. Thus, the 

20 motion state of the vehicle may not be accurately represented 
in such a neural network. In particular, the values output 
from the neural network (for example, the side slipping angle 
and the yaw rate) are varied in accordance with not only the 
input, but also the value thereof at the present time (a present 

25 value) . Consequently, it is necessary to feed back the output 
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value and reflect the output value to the neural network, in 
order to estimate the motion state of the vehicle with high 
precision. In this sense, the JP-A-6-286630 supplies the time 
delay value of the output value to an input layer by using an 
5 ARAM model, thereby enhancing the estimation precision of the 
road surface friction coefficient . However , the neural network 
having such feedback has a problem that the coupling weight 
coefficient cannot be learned according to the principle of 
a learning rule such as back propagation. Thus, accurate 
10 estimation of the road surface friction coefficient is hardly 
achieved. 

The present invention has been implemented in view of 
the foregoing situation, and has an object to provide a novel 
method of creating amotion model of a vehicle by using a recurrent 
15 neural network containing a feedback loop. 

SUMMARY OF THE INVENTION 

The present invention has been implemented in view of 
the foregoing situation, and has an object to provide a novel 
20 method of creating amotion model of a vehicle by using a recurrent 
neural network containing a feedback loop. 

In order to solve such problems , a first invention provides 
a vehicle motion model generating device for generating a vehicle 
motion model that represents a motion state of a vehicle. This 
25 vehicle motion model generating device includes a first 
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recurrent neural network formed by connecting plural nodes so 
that an output of a node is input to another node in accordance 
with a predetermined coupling weight coefficient, including 
a loop feeding back an output of at least one node to at least 
5 one of the one node and a node other than the one node, and 
an optimizing unit for determining an optimum solution of the 
coupling weight coefficient in the first recurrent neural 
network based on a learning rule using a hereditary algorithm. 
Here, the first recurrent neural network outputs a first 

10 parameter indicating the motion state of the vehicle based on 
predetermined input information, thereby functioning as the 
vehicle motion model. 

In the first invention , the first recurrent neural network 
may have a hierarchical structure including at least an input 

15 layer formed of one or more nodes and an output layer formed 
of one or more nodes. In this case, it is preferable that the 
optimizing unit determines the optimum solution of the coupling 
weight coefficient with connection of respective nodes between 
neighboring layers being set as a processing target. 

20 Correspondingly, in the first invention, the first recurrent 
neural network may be formed of the plural nodes connected 
mutually so that the output of the one node is input to all 
the plural nodes including the one node, and respective outputs 
of the plural nodes are input to the one node. In this case, 

25 it is preferable that the optimizing unit determines the optimum 
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solution of the coupling weight coefficient with mutual 
connection of the plural nodes being set as a processing target . 

In the first invention, it is preferable that each of 
the plural nodes uses one of a sigmoid function and a non- sigmoid 
5 function other than the sigmoid function as a transfer function . 

Further, in the first invention, the vehicle motion state 
model generating device may includes a second recurrent neural 
network constructed as a network different from the first 
recurrent neural network, functioning as the vehicle motion 

10 model by outputting a second parameter indicating a motion state 
of the vehicle different from the first parameter. It is 
preferable that the optimizing unit further determines the 
optimum solution of the coupling weight coefficient in the second 
neural network based on the learning rule using the hereditary 

15 algorithm. In this case, it is preferable that the first 
recurrent neural network and the second recurrent neural network 
are mutually connected to each other so that a state variable 
having a correlation with the first parameter output from the 
first recurrent neural network is input to the second neural 

20 network. Here, the state variable represents a road surface 
state or a motion state of the vehicle. 

A second invention provides a method for generating a 
vehicle motion model that represents a motion state of a vehicle . 
The method, in which a first recurrent neural network is formed 

25 by connecting plural nodes so that an output of a node is input 
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to another node in accordance with a predetermined coupling 
weight coefficient and includes a loop feeding back an output 
of at least one node to at least one of the one node and a node 
other than the one node, is executed by a computer, and includes 
5 the steps of determining an optimum solution of a genetic type 
based on a learning rule using a hereditary algorithm while 
setting the coupling weight coefficient in the first recurrent 
neural network as the genetic type, and outputting an optimum 
solution of the coupling weight coefficient to the first 

10 recurrent neural network based on the optimum solution of the 
genetic type . Here , the first recurrent neural network outputs 
a parameter indicating the motion state of the vehicle based 
on predetermined input information, thereby functioning as the 
vehicle motion model. 

15 In the second invention, the first recurrent neural 

network may have a hierarchical structure including at least 
an input layer formed of one or more nodes and an output layer 
formed of one or more nodes. Here, the method may include the 
step of determining the optimum solution of the genetic type 

20 while connection of respective nodes between neighboring layers 
is set as a processing target. Correspondingly, the first 
recurrent neural network may be formed of the plural nodes 
connected mutually so that the output of the one node is input 
to all the plural nodes including the one node, and respective 

25 outputs of the plural nodes are input to the one node. Here, 
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the method may include the step of determining the optimum 
solution of the coupling weight coefficient while mutual 
connection of respective nodes is set as a processing target. 

In the second invention, it is preferable that a second 
recurrent neural network is constructed as a network different 
from the first recurrent neural network, and functions as the 
vehicle motion model by outputt ing a second parameter indicating 
amotion state of the vehicle different from the first parameter . 
The method preferably includes the steps of determining the 
optimum solution of the genetic type while setting the coupling 
weight coefficient in the second recurrent neural network as 
the genetic type, and outputting the optimum solution of the 
coupling weight coefficient to the second recurrent neural 
network based on the optimum solution of the genetic type. 
Further , it is preferable that the first recurrent neural network 
and the second recurrent neural network are mutually connected 
to each other so that a state variable having a correlation 
with the first parameter output from the first recurrent neural 
network is input to the second neural network. Here, the state 
variable represents one of a road surface state and a motion 
state of the vehicle. 

A third invention provides a road surface friction 
coefficient estimating device for estimating a road surface 
friction coefficient. This road surface friction coefficient 
estimating device estimates the road surface friction 
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coefficient based on a vehicle motion model that represents 
a motion state of a vehicle and is generated by a vehicle motion 
model generating device which is provided in the above first 
invention . 

5 A fourth invention provides a vehicle behavior estimating 

device for estimating a behavior of a vehicle. This vehicle 
behavior estimating device estimates the behavior of the vehicle 
based on a vehicle motion model that represents a motion state 
of the vehicle and is generated by a vehicle motion model 
10 generating device that is provided in the above first invention . 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing the construction of 
a vehicle motion model generating device according to an 
15 embodiment ; 

Fig. 2 is a diagram showing the structure of RNN which 
functions as an estimating module; 

Fig. 3 is a flowchart showing the procedure of determining 
the optimum solution of a weight coefficient and a threshold 
20 value by using a genetic algorithm; 

Fig. 4 is a diagram showing the concept of an estimating 
function; 

Fig. 5 is a diagram showing a learning result concerning 
a first -generation individual of the RNN; 
25 Fig. 6 is a diagram showing a learning result concerning 
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a 100-th generation individual of the RNN; 

Fig. 7 is a diagram showing a learning result concerning 
a 30, 000-th generation individual of the RNN; 

Fig, 8 is a diagram showing an example of a condition 
5 given to the RNN; 

Fig. 9 is a diagram showing the network structure of the 
RNN constituting each estimating module according to a second 
embodiment ; 

Fig. 10 is a diagram showing the corresponding 
10 relationship in coupling weight coefficient among nodes; and 
Fig. 11 is a diagram showing an example of a transfer 
function. 

DESCRIPTION OF PREFERRED EMBODIMENTS 

15 (First Embodiment) 

Fig. 1 is a block diagram showing the construction of 
a vehicle motion model generating device for generating a vehicle 
motion model according to the present embodiment . This vehicle 
motion model generating device or, simply a generating device 

20 generates a vehicle motion model representing the motion state 
of a vehicle or a vehicle motion state. The generating device 
1 may be used as a computer comprising a CPU, a RAM, a ROM, 
an input /output interface, etc. As not clearly shown in Fig. 
1, an input device containing a keyboard and a mouse, a display 

25 device such as CRT and a liquid crystal display device are 
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connected to the generating device 1 . An operator can input 
necessary information, etc. by operating the input device based 
on the information displayed on the display device. From the 
functional aspect of the generating device 1, it comprises an 
5 optimizing unit 10 and a vehicle motion model unit 20. 

The optimizing unit 10 carries out an operation based 
on a learning rule using a genetic algorithm, and outputs 
predetermined parameters as an operation result to the vehicle 
motion model unit 20. As described later, these parameters 

10 correspond to a coupling weight coefficient Ki j of a recurrent 
neural network and a threshold value 6 j , since each of estimating 
modules 21 to 24 constituting the vehicle motion model unit 
20 is constructed by the recurrent neural network. Here, "the 
recurrent neural network" means a neural network having a 

15 feedback loop. 

The vehicle motion model unit 20 includes one or more 
estimating modules (in the present embodiment, the estimating 
modules 21 to 24), where each of the estimating modules 21 to 
24 functions as a vehicle motion model representing a distinct 

20 motion state. Each of the estimating modules 21 to 24 includes 
a recurrent neural network (hereinafter merely referred to as 
" RNN " ) , and outputs a parameter indicating the motion state 
of a vehicle (hereinaf tier referred to as "vehicle parameters" ) 
based on predetermined input information. In the present 

25 embodiment , four vehicle parameters are output from the vehicle 
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motion model unit 20. Specifically, the estimating module 21 
outputs an estimation value of the yaw rate as a vehicle parameter , 
and the estimating module 22 outputs an estimation value of 
the lateral acceleration (hereinafter merely referred to as 
5 "lateral G") as a vehicle parameter. Furthermore, the 
estimating module 23 outputs an estimation value of roll as 
a vehicle parameter, and the estimating module 24 outputs an 
estimation value of pitch as a vehicle parameter. 

Fig. 2 is a diagram showing the structure of the RNN which 

10 functions as the estimating module . Here, the RNN constituting 
the estimating module 21 for the yaw rate is used for description 
of the RNN. On the other hand, the description of the other 
estimation modules 22 to 24 is omitted because the RNNs 
constituting the other estimating module 22 to 24 have basic 

15 structures similar to that of the estimating module 21. In 
the RNN, a network is formed by connecting plural nodes Nn For 
example, nodes Nl to N18 are shown in Fig. 2. The plural nodes 
Nn are hierarchically sectioned, thereby forming a hierarchical 
structure including an input layer, an intermediate layer and 

20 an output layer. The number of the nodes Nn constituting each 
layer is properly set by the operator . In this case , the operator 
determines the structure of the network so that each layer has 
the optimum node number in consideration of the conflicting 
relationship between the enhancement in reliability of the 

25 solution achieved by increasing the number of the nodes Nn and 
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enhancement of processing speed achieved by reducing the number 
of the nodes Nn . 

The output layer is equipped with one node Nn ( the node 
N16) for output ting a vehicle parameter. However, in the 
5 present embodiment, a feedback loop from the output layer to 
the input layer exists as described later. Additionally, the 
node N17 aiming only at the feedback is further set to carry 
out the feedback to the input layer in a multiplexing style. 
Furthermore, in Fig. 2, the node N18 for outputting a state 

10 variable P is set separately from the vehicle parameter in the 
RNN constituting the estimating module 21 for the yaw rate. 
The state variable P output from the node N18 is input to the 
input layer node Nn of the RNN constituting the estimating 
modules 22 to 24 for the lateral G, the roll and the pitch. 

15 The numbers of nodes Nn necessary to input the input 

information are set in the input layer. Here, the input 
information is required for the RNN to output the vehicle 
parameter. In the present embodiment, the input information 
corresponds to a steering angle, a steering angular velocity, 

20 a steering angular acceleration, a steering reaction force, 
a vehicle speed, a vehicle acceleration, etc. Furthermore, 
the nodes N7 to N10 associated with the feedback are also set 
in the input layer, in connection with the feedback from the 
outputs of the output layer nodes N16, N17 to the input layer. 

2 5 The nodes Nn for inputting the state variable P shown in Fig. 
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1 are further set in the other estimating modules 22 to 24. 
These nodes Nn are not specifically shown in Fig. 2 because 
the RNN constituting the estimating module 21 for the yaw rate 
is illustrated. 

Considering the above-described relationship , a required 
number of the nodes Nn are prepared based on an operator's 
experience, for instance. In the example shown in Fig. 2, the 
intermediate layer includes five nodes Nn of the nodes Nil to 
N15. 

When the nodes Nn constituting each layer is set as 
described above , the connection among the nodes is automatically 
determined. In the hierarchical type RNN, the connection of 
the among is determined so that the output of each node Nn 
constituting some layer is input to all the nodes Nn constituting 
the layer at the subsequent stage. For example, the node Nl 
of the input layer is input to all the nodes Nil to N15 of the 
intermediate layer as shown in Fig. 2. 

As one feature of the present embodiment, the RNN has 
a loop for feeding back the output of some node Nn . Specifically , 
as shown in Fig. 2, the output of the output layer node N16 
is output as a yaw rate which is one of the vehicle parameters , 
delayed for only one cycle by a first delay circuit Z" 1 , and 
then, input to the input layer node N7 . The output from the 
first delay circuit Z" 1 is further delayed for only one cycle 
by a second delay circuit Z" 1 equipped at the subsequent stage 
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thereto (accordingly, after delayed for two cycles with respect 
to the initial input corresponding to the output of the node 
N16) , and then input to the input layer node N8 . On the other 
hand, the output of the node N17 of the output layer is delayed 
5 for only one cycle by a third delay circuit Z" 1 different from 
the first and second delay circuits Z" 1 , and then input to the 
input layer node N9 . Furthermore, the output from the third 
delay circuit Z" 1 is further delayed for only one cycle by a 
fourth delay circuit Z" 1 equipped at the subsequent stage 

10 (accordingly, after delayed for only two delays with respect 
to the initial input corresponding to the output of the node 
N17), and then input to the input layer node N10. 

In such a RNN, each node Nn generates one output Yi for 
a predetermined number of inputs yi according to a predetermined 

15 transfer function. The output Yi is output to the nodes Nn 
constituting the layer at the subsequent stage. Specifically, 
each node Nn carries out the calculation represented equations 
1, 2 on the input data yi, and outputs the calculation result 
as an output data Yj . Here, Kij represents the coupling weight 

20 coefficient between an (n=i)-th (i=l to 15) node Ni and an 
(n=j)-th (j = ll to 15 (i = l to 10), j = 16 to 18(i=ll to 15)) node 

N j , and 9j represents a threshold value. 
[Equation 1] 

Xj = 2Kij • yi(i = 1 to n) 
25 [Equation 2] 



15 



Yj = 1/(1 + exp (- (Xj - 83))) 

As apparent from the above equations 1, 2, the output 
of some node Nn is input to another node Nn in accordance with 
a predetermined coupling weight coefficient Ki j . Here, the 
5 equation 2 is called as a sigmoid function which generally used 
as a function for nodes in the RNN . The sigmoid function 
continuously varies from 0 to 1, and the threshold value 0j 
approaches to a step function as it is reduced. 

When the vehicle parameter is estimated by using the RNN, 

10 proper adjustment (learning) of the coupling weight coefficient 
Ki j and the threshold value 9 j is required in order to enhance 
the precision of the estimation result. The optimum solution 
of the coupling weight coefficient Ki j and the threshold value 
0j is learned and properly determined by the optimizing unit 

15 10 described above. The connection to be processed by the 
optimizing unit 10 becomes the connection of the respective 
nodes between the neighboring layers except for the feedback 
loop. 

Fig. 3 is a flowchart showing the procedure of determining 
20 the optimum solution of the weight coefficient wij and the 
threshold value 0j by using the genetic algorithm. The 
processes shown in Fig. 3 are carried out by the optimizing 
unit 10. The optimizing unit 10 carries out the following 
processes on each of the estimation modules 21 to 24. In this 
25 case, the estimation module (RNN) 21 for the yaw rate will be 
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described as an example. Since the same processes are carried 
out on the other estimation modules (RNN) 22 to 24, the 
description of the other estimation modules 22 to 24 is omitted. 

Initially, in step SI, an initial individual assembly 
5 is formed by a set of N individuals each including a genetic 
type . Here , the genetic type corresponds to the coupling weight 
coefficient Kij and the threshold value 6 j . In this case, one 
individual includes a genetic type containing all the coupling 
weight coefficients Kij and threshold values 0j in the RNN. 

10 The initial individual assembly is a set of N individuals 
described above, and includes genetic types Kij, 8j varied by 
each individual. For example, the initial value of the genetic 
type Kij , 9 j constituting each individual 1 to N is determined 
by a random number. In other words, if a coupling weight 

15 coefficient Kij and a threshold value 9j are substituted into 
the RNN for every individual 1 to N, the RNNs of N are calculated, 
each having the determined coupling weight coefficients Kij 

and the threshold values 0 j . 

The next, in step S2 , an adaptive value A is calculated. 

20 On the assumption that the adaptive value A is calculated, the 
operator beforehand obtains through experiments or simulations 
a vehicle parameter which is actually indicated by a vehicle 
under a predetermined condition (hereinafter referred to as 
"an actual vehicle parameter"). The operator stores a map 

2 5 describing the relationship between the predetermined 
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condition and the actual vehicle parameter into the ROM of the 
generating device 1. The predetermined condition described 
in the map is parameters input to the input layer node Nn except 
for the feedback elements , and corresponds to the steering angle , 
5 steering angular velocity, the steering angular acceleration, 
the steering reaction force, the vehicle speed and the vehicle 
acceleration. Based on the above-described assumption, 
referring to the map, the optimizing unit 10 inputs these 
conditions to the respective input layer nodes Nl to N6 of the 

10 RNN, and acquires an estimation value of the vehicle parameter 
(hereinafter referred to as "an estimated vehicle parameter") 
as the output of the output layer node N16. The adaptive value 
A is calculated based on the estimated vehicle parameter and 
the actual vehicle parameter by using an estimation function 

15 represented by an equation 3. 
[Equation 3] 
A = 2|Vt - Vc| 

Fig. 4 is a schematic diagram showing the estimation 
function. In Fig. 4, time- sequential variation of an actual 

20 vehicle parameter (for example, an actual yaw rate ) is indicated 
by a solid line, and time- sequential variation of an estimated 
vehicle parameter (for example, an estimated yaw rate) is 
indicated by a broken line. As shown in the equation 3, the 
estimation function corresponds to the integration of 

25 difference (an absolute value) between the actual vehicle 
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parameter Vt and the estimated vehicle parameter Vc (the area 
of a hatched area shown in Fig. 4) . In other words , the estimation 
function is set so that the adaptive value A is lowered, as 
the estimated yaw rate is closer to the actual yaw rate. 
5 The optimizing unit 10 calculates the adaptive value A 

for each individual 1 to N, and rearranges the individuals 1 
to N in the increasing order of the adaptive value A. As a 
result, the individuals 1 to N of the individual assembly are 
rearranged as individuals 1 to N, in order of increasing 

10 difference between the actual vehicle parameter and the 
estimated vehicle parameter. 

In step S3 f it is judged whether the value of the adaptive 
value A concerning the individual 1 is greater than the value 
of a judgment adaptive value Aerror. The judgment adaptive 

15 value Aerror is set as the maximum value of the adaptive value 
A at which the estimated yaw rate and the actual yaw rate can 
be regarded as being substantially identical to each other. 
Accordingly, if a positive judgment is made in the step S3 (the 
adaptive value A is greater than the judgment adaptive value 

20 Aerror), the process shifts to the subsequent step S4 . Else, 
if a negative judgment is made in the step S3 (the adaptive 
value A is equal to or smaller than the judgment adaptive value 
Aerror) , the process shifts to step S6 at the subsequent stage. 

In the step S4, the individuals 1 to N are selected or 

2 5 eliminated by use of a known manner. Such a 
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selecting/eliminating manner is, for instance, inverse 
roulette type selection, ranking type selection, tournament 
type selection, or the like. Through the process of the step 
S4 , an individual or an individual group is selected, and deleted 
5 from the individual assembly. The same numbers of individuals 
are shifted to the position at which individuals are deleted 
as described above, for example, in the increasing order of 
the adaptive value A, whereby the number of the individuals 
constituting the population is maintained. 

10 In step S5, a next generation population is generated. 

Specifically, the optimizing unit 10 carries mutation and 
chiasma to the genetic type Kij , 9 j constituting an individual 
in the assembly. The mutation means that any genetic type Kij 
(or 0 j ) of some individual is selected by a random number or 

15 the like, then, converted to a value generated by another random 
number. Furthermore, the chiasma means that the values of 
genetic types Kij (or0j) selected by random numbers are mutually 
exchanged to each other, in some individual group selected by 
a random number. However, the method of selecting an 

20 individual ( s ) subjected to the mutation and the chiasma is not 
limited to the above-described random number selecting method. 
Besides , with respect to an individual having the lowest adaptive 
value A, it may not be subjected to mutation and chiasma to 
keep the individual concerned. The process returns to the step 

25 S2 and the subsequent steps to repeat the above-described 
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processes until the adaptive value A of the individual 1 becomes 
equal to or lower than the adaptive value Aerror. 

On the other hand, in the step S6 associated with the 
positive judgment of the step S3, the genetic types Ki j , 0j 
5 constituting the individual 1 are determined as the optimum 
solution, and the process exits the entire routine. In this 
case , the optimizing unit 10 determines the values of the genetic 
types Ki j , 9 j as the optimum solution of the corresponding 
coupling weight coefficients Ki j and threshold values 6j based 

10 on the genetic types Kij , 6j constituting the individual 1. 
The optimum solution thus determined is output to the RNN, and 
set to the coupled coefficient Kij and the threshold value 6j 
of the RNN in each of the estimating modules 21 to 24. 
Consequently, the vehicle motion model unit 20 containing the 

15 estimation modules 21 to 24 alone or the generating device 1 
itself functions as a vehicle motion model. 

Figs. 5 to 7 show learning results of the RNN. Fig. 5 
shows a learning result concerning a first-generation 
individual. Fig. 6 shows a learning result concerning 100-th 

20 individual. Fig. 7 shows alearning result concerning 30 , 000- th 
individual. In these figures, an actual vehicle parameter (a 
yaw rate) under the above -described predetermined condition 
(for example, the steering angle shown in Fig. 8) is drawn as 
a thin solid line. Meanwhile, the estimated vehicle parameter 

2 5 (the yaw rate) under the same condition is drawn as a heavier 
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solid line. As apparent from these figures, in the first 
generation , both the lines show the same tendency in an amplitude 
period, however, are greatly different in value. Thus, the 
estimated vehicle parameter does not reproduce the actual 
5 vehicle parameter. In the 100-th generation , although both 
the lines show the same tendency, a difference between the lines 
appears in the neighborhood of the local maximum (or minimum) 
point. In the 30, 000-th generation, both the lines show 
substantially the same value over the entire area, where the 

10 thin line corresponding to the actual vehicle parameter is 
covered by the heavy solid line corresponding to the estimated 
vehicle parameter. As described above, by learning until the 
30,000-th generation or the like, the estimated parameter can 
reproduce the actual vehicle parameter with high level or 

15 precision. That is, by learning with the genetic algorithm, 
the optimum solution of the weight coefficient Kij of the RNN 
containing the feedback loop can be effectively derived. 

As described above, the learning of the coupling weight 
coefficient Kij (further, the threshold value 6 j ) of the RNN 

20 is carried out by using the genetic algorithm to determine the 
optimum solution, according to the present embodiment. The 
RNN containing feedback has a problem not able to use a learning 
rule based on the principle of a steepest lowering method such 
as back propagation. However, this problem can be solved by 

25 using the genetic algorithm, according to the present 
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embodiment . 

By using the RNN , the output of the output layer node 
Nn containing the vehicle parameter is fed back to the input 
layer node Nn. Therefore, the value fed back to the RNN is 
5 reflected, so that the RNN of the present embodiment can estimate 
the vehicle parameter while the previous output thereof is added 
to the estimation. As a result, the method according to the 
present embodiment can generate a vehicle motion model which 
represents the motion state of an actual vehicle more accurately 

10 than the f eed- forward type neural network, that is, has excellent 
reproducibility of the motion state of a vehicle. 

In the present embodiment , each of the estimating modules 
21 to 24 outputs one vehicle parameter. For example, a node 
Nn for outputting the above four vehicle parameters may be 

15 equipped to the output layer in one RNN, whereby the one RNN 
can be substituted for four RNNs. However, according to the 
system construction of the present embodiment, the number of 
the nodes Nn constituting each RNN can be reduced by designing 
a node Nn for outputting a vehicle parameter as a single body. 

20 Accordingly, the learning speed of the optimizing unit 10 is 
increased for each of the estimating modules 21 to 24. 
Furthermore, as the learning speed is increased, the number 
of individuals in an initial individual assembly could be 
increased. Thus, learning can be performed while various 

2 5 individuals are contained in the individual assembly when the 
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learning is carried out by using a genetic algorithm. 
Accordingly, the optimum solution of the genetic type Ki j , 0j 
can be determined with high precision. Additionally, if the 
number of nodes is small, there could be achieved such an effect 
5 that the combination of the nodes can be established in a 
small-scale, so that the number of dimension in solution space 
in the genetic algorithm can be reduced. In conclusion, the 
optimum solution is more liable to be detected. 

The vehicle parameters output from the estimating modules 

10 21 to 24 represent the motion state concerning the common vehicle . 
Therefore, even if the types thereof are different, the vehicle 
parameters must vary synchronously to each other due to the 
motion of the vehicle. Therefore, according to the present 
embodiment, the state variable P is output from the estimating 

15 module 21 for the yaw rate, and then, input to the estimating 
modules 22 to 24, for instance. The state variable P output 
from the estimating module 21 for the yaw rate has high 
correlation with the vehicle parameter (specifically, the yaw 
rate), that is, the output result of the estimating module 21 

20 by learning the coupling weight coefficient Kij with use of 
the genetic algorithm. For example, the state variable P 
represents the motion state of the vehicle or the road surface 
state corresponding to the yaw rate of the vehicle. Therefore, 
in the estimating modules 2 2 to 24 to which the state variable 

25 P is input , the coupling weight coefficient Kij is learned under 



the condition that the state variable P having the correlation 
with the yaw rate is included, so that the RNNs constituting 
the four estimating modules 21 to 24 can be synchronized one 
another. Basically, the state variable P functions as a 
5 reference value for synchronization with the estimating module 
21 in the estimating modules 22 to 24. Accordingly, in the 
system construction of the present embodiment, each vehicle 
parameter can reproduce a value closer to the motion state of 
the actual vehicle as compared with a case where the estimating 

10 modules 21 to 24 are constructed independently of one another. 

The state variable P functioning as described above is 
not limited to the output from the estimating module 21 for 
the yaw rate, and may be any one of the outputs of the estimating 
modules 22 to 24. However, among the vehicle parameters, the 

15 yaw rate can be detected with relatively higher precision than 
the other vehicle parameters. Thus, the state variable P is 
preferably output from the estimating module 21 for the yaw 
rate (in short, the reliability of the state variable P is 
enhanced) . Furthermore, the present embodiment may be modified 

20 so that each state variable P is output from the estimating 
modules 21 to 24, then, input to an estimating module different 
from an estimating module at the output side, thereby making 
the estimating modules have association ( synchronism) with each 
other. In the above- described embodiment, the state variable 

2 5 P is handled as a variable which has correlation with the yaw 
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rate and represents the road surface state or the motion state 
of the vehicle. However, from the viewpoint that the state 
variable P has correlation with the yaw rate, the vehicle 
parameter (the yaw rate) output from the node N16 itself may 
5 be set as the state variable P. 

In the present embodiment, the vehicle motion model unit 
20 is constructed by the four estimating modules 21 to 24. 
However, the invention is not limited to this mode. The number 
of estimating modules may be arbitrarily determined in 

10 accordance with the necessary vehicle parameters . In addition , 
the vehicle parameters are not limited to the four values, which 
are, the yaw rate, the lateral G, the roll and the pitch, and 
may include other parameters such as a vehicle slipping angle, 
a front -wheel slipping angle, a rear-wheel slipping angle, the 

15 friction coefficient between the wheels and the road surface, 
etc . 

Furthermore, in the present embodiment, the coupling 
weight coefficient Kij and the threshold value 6j of the RNN 
are learned by referring to a map and comparing this map with 

20 the output of the RNN. However, the generating device 1 may 
be mounted in an actual vehicle to carry out the learning on 
a real-time basis. In this case, the steering angle, the 
steering angular velocity, the steering angular acceleration, 
the steering reaction force, the vehicle speed, and the vehicle 

25 acceleration which correspond to the inputs are detected from 



sensors or calculated based on detected values, and then, set 
as inputs of the RNN . The vehicle parameters corresponding 
to the outputs are detected by sensors, and the detected values 
and the outputs of the RNN are compared. 
5 (Second Embodiment) 

Fig. 9 is a diagram showing the network structure of the 
RNN constituting each of the estimating modules 21 to 24 
according to a second embodiment. In the embodiment shown in 
Fig. 9, some nodes Nn are omitted from the illustration for 

10 the same of convenience. Each of the estimating modules 21 
to 24 corresponds to the RNN comprising the plural nodes Nn 
and the optimum solution of the coupling weight coefficient 
Ki j and the threshold value 0 j is determined by using the genetic 
algorithm, similarly to the first embodiment. However, 

15 according to the second embodiment, the network structure of 
the RNN serving as a premise is different from that of the first 
embodiment. The RNN according to the first embodiment has a 
hierarchical structure, where the connection among the nodes 
Nn is determined in advance by setting the number of the nodes 

20 Nn constituting each layer. On the other hand, in the RNN 
according to the second embodiment , plural nodes are mutually 
connected to each other so that some node Nn is coupled to all 
nodes including the node Nn itself. Accordingly, the output 
of the node Nl is input to all the nodes Nl to Nn including 

25 the node Nl itself. Meanwhile, the outputs of all the nodes 
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Nl to Nn are input to the node Nl. 

Fig . 10 is a diagram showing the correspondence of coupling 
weight coefficient Kij to each other among the nodes. In the 
network as described above, the coupling weight coefficients 
Kij among the nodes are represented as a matrix shown in Fig. 
10. This matrix represents the coupling weight coefficients 
Kij concerning the connection among all the nodes irrespective 
of the type such as feedforward and feedback. In the present 
embodiment, the optimizing unit 10 in the above-described 
network structure determines the optimum solution of the 
coupling weight coefficient Kij and the threshold value 0 j , 
with the mutual connection of plural nodes as a processing 
target . 

The method of determining the optimum solution of the 
coupling weight coefficient Kij and the threshold value 9 j is 
basically the same as in the first embodiment. Therefore, only 
a brief description thereof will be made hereunder, and the 
detailed description is omitted. Specifically, first, the 
coupled coefficient Kij (that is, Kll to Knn) and the threshold 
value 6j shown in Fig. 10 are set to a genetic type, and plural 
individuals comprising genetic types having various values are 
prepared (an initial individual assembly is generated). The 
adaptive value A is calculated for every individual. In this 
case, random six nodes Nn are selected from the nodes Nn 
constituting the network. In addition, the steering angle, 
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the steering angular velocity, the steering angular 
acceleration, the steering reaction force, the vehicle speed 
and the vehicle acceleration are input to the selected nodes 
based on the above -described map. The next, any one of the 
5 nodes Nn different from the nodes used for input is selected, 
and the output is acquired from the node Nn thus selected . Based 
on this output and the vehicle parameters described on the map, 
the adaptive value A is calculated by using the estimation 
function shown by the equation 3. Similarly to the processes 
10 shown in Fig. 3, the selection, elimination and generation of 
the next generation population are repeated until the adaptive 
value A becomes equal to or smaller than the judgment adaptive 

value Aerror. Consequently, the optimum genetic type Kij, 8j 
is determined, whereby the optimum solution of the coupling 
15 weight coefficient Kij and threshold value 0j in the RNN is 
determined. 

As described above, the present embodiment adopts the 
RNN having the network structure whose coupling state is not 
initially limited or fixed , and in which feedforward and feedback 

20 are contained and all the nodes are connected to each other. 
Even in such a network structure , the coupling weight coefficient 
Kij (and threshold value 0 j ) can be determined by using the 
genetic algorithm. In this case, if the value of the determined 
coupling weight coefficient Kij becomes "0", coupling of the 

25 node Ni and the node Nj becomes redundant. In other words. 



29 



although all the nodes Nn are formally coupled one another as 
a network, only necessary nodes remain substantially in 
accordance with the value of the coupling weight coefficient 
Kij. More specifically, by learning the coupling weight 
5 coefficient Kij according to the genetic algorithm, nodes 
necessary to calculate the vehicle parameter are equivalently 
selected. Of course, it is expected as a learning result that 
the network structure as shown in Fig. 2 be generated in 
accordance with the value of the coupling weight coefficient 

10 Kij. However, in such a case, only the necessary nodes are 
selected from the RNN in which all the nodes Nn are connected 
to each other. Thus, it is advantageously guaranteed by the 
genetic algorithm that the network structure concerned is the 
optimum structure to estimate the vehicle parameter. 

15 In the first and second embodiments, the foregoing 

description has been made on the assumption that the transfer 
function is the sigmoid function. However, if the learning 
is carried out by using the genetic algorithm, the transfer 
function of the nodes constituting the RNN is not limited to 

20 the sigmoid function. For example, a non- sigmoid function as 
shown in Fig. 11 may be used as the transfer function. Fig. 
11(a) shows an example of a function in which the output indicates 
a non-linear characteristic called as hysteresis with respect 
to the input. Figs. 11(b) to (d) show examples of functions 

25 indicating non-linear characteristics called binarization , 
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dead zone and limiter. 

The transfer function is not limited to the functions 
indicating these non-linear characteristics, and may be a 
function containing an element like a sequence circuit. That 
5 is , a function in which the output is determined by a past input 
sequence may be used as the transfer function . As such a function 
may be used time delay (Z" 1 ) , peak hold (the peak value during 
a predetermined interval is output), maximum value (the 
time-sequential peak value is output), integration, 

10 differentiation, etc. 

By including non- sigmoid functions in the transfer 
function of the nodes Nn as described above, a vehicle motion 
model that can represent even input/output characteristics 
which cannot be represented by the sigmoid function can be 

15 generated. Accordingly, the vehicle motion model as described 
above can represent an input /output response when a rear wheel 
of a vehicle slips , thus , the motion model of the vehicle changes 
drastically at a dash. Therefore, the reproducibility of the 
motion state of the vehicle can be enhanced. 

20 The vehicle motion model generating device for generating 

the vehicle motion model as described above or the vehicle motion 
model generated by the vehicle motion model generating device 
is usable while installed in a part of a road surface friction 
coefficient estimating device for estimating the road surface 

25 friction coefficient or a vehicle behavior estimating device 
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for estimating the behavior of a vehicle. For example, when 
this vehicle motion model is applied to the road surface friction 
coefficient estimating device, there are prepared a high-ji 
vehicle motion model in which the motion state of the vehicle 
5 on the road surface having a high road surface friction 
coefficient has been learned. Additionally, a low-^i vehicle 
motion model in which the motion state of the vehicle on the 
road surface having a low road surface friction coefficient 
has been learned. Thereafter, the steering angle, the steering 

10 angular velocity, the steering angular acceleration, the 
steering reaction force, the vehicle speed, and the vehicle 
acceleration concerning the present vehicle are input to each 
vehicle motion model to derive a vehicle parameter such as the 
yaw rate. Furthermore, the motion state of the vehicle 

15 corresponding to the vehicle parameter thus derived is detected 
by a sensor or the like. This detected value is compared with 
the value derived from the vehicle motion model. Accordingly, 
the road surface state corresponding to the learning condition 
of one vehicle motion model exhibiting such a tendency that 

20 the above values are approximated to each other can be estimated 
as the present road surface friction coefficient . On the other 
hand, the vehicle behavior estimating device inputs the steering 
angle, the steering angular velocity, the steering angular 
acceleration, the steering reaction force, the vehicle speed 

25 and the vehicle acceleration concerning the present vehicle 



to the vehicle motion model to derive the vehicle parameter 
such as the yaw rate, and estimates a behavior the vehicle 
indicates. Furthermore, by using this vehicle motion model, 
the behavior of the vehicle and ride quality of the vehicle 
5 can be also estimated. 

As described above, according to this invention, the 
optimum solution of the coupling weight coefficient of the 
recurrent neural network having the feedback loop can be 
determined by using the genetic algorithm. Furthermore, the 

10 feedback element can be taken into consideration in the recurrent 
neural network by applying the recurrent neural network to the 
vehicle motion model. Accordingly, there can be generated a 
vehicle motion model more suitable for the motion state of the 
actual vehicle than the feed- forward type recurrent neural 

15 network, or enhanced in reproducibility of the motion state 
of the vehicle. 

While the present invention has been disclosed in terms 
of the preferred embodiments in order to facilitate better 
understanding of the invention, it should be appreciated that 

20 the invention can be embodied in various ways without departing 
from the principle of the invention. Therefore, the invention 
should be understood to include all possible embodiments which 
can be embodied without departing from the principle of the 
invention set out in the appended claims. 

25 Additionally, the disclosure of Japanese Patent 



Application No. 2003-041507 filed on February 19 , 2003 including 
the specification, drawing and abstract is incorporated herein 
by reference in its entirety. 
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